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The present-day Universe is highly magnetized, even though the first magnetic seed fields were 
most probably extremely weak. To explain the growth of the magnetic field strength over many 
orders of magnitude fast amplification processes need to operate. The most efficient mechanism 
known today is the small-scale dynamo, which converts turbulent kinetic energy into magnetic energy 
leading to an exponential growth of the magnetic field. The efficiency of the dynamo depends on 
the type of turbulence indicated by the slope of the turbulence spectrum v(£) oc , where v(£) is the 
eddy velocity at a scale I. We explore turbulent spectra ranging from incompressible Kolmogorov 
turbulence with # = 1/3 to highly compressible Burgers turbulence with ■& = 1/2. In this work 
we analyze the properties of the small-scale dynamo for low magnetic Prandtl numbers Pm, which 
denotes the ratio of the magnetic Reynolds number, Rm, to the hydrodynamical one, Re. We solve 
the Kazantsev equation, which describes the evolution of the small-scale magnetic field, using the 
WKB approximation. In the limit of low magnetic Prandtl numbers the growth rate is proportional 
to Rni''"''"'^'. We furthermore discuss the critical magnetic Reynolds number Rm cr i t , which 
is required for small-scale dynamo action. The value of Rm cr it is roughly 100 for Kolmogorov 
turbulence and 2700 for Burgers. Furthermore, we discuss that Rm cr i t provides a stronger constraint 
in the limit of low Pm than it does for large Pm. We conclude that the small-scale dynamo can 
operate in the regime of low magnetic Prandtl numbers, if the magnetic Reynolds number is large 
enough. Thus, the magnetic field amplification on small scales can take place in a broad range of 
physical environments and amplify week magnetic seed fields on short timescales. 



I. INTRODUCTION 

A large fraction of the Universe is magnetized. Various 
astrophysical phenomena have their origin in strong mag- 
netic fields, for example jets from stars or galaxies and 
stellar activity. The question arises where those strong 
fields came from, especially because the generation mech- 
anisms during inflation pQ , phase transitions in the early 
Universe [2] or battery processes [3H5] typically produce 
very week seed fields. 

Magnetohydrodynamical dynamos are the most efficient 
mechanisms known to amplify weak magnetic seed fields. 
In particular the small-scale or turbulent dynamo is im- 
portant as it converts turbulent kinetic energy into mag- 
netic energy on very short timescales. 
After the magnetic seed fields have been amplified expo- 
nentially by the kinematic dynamo on small scales, the 
nonlinear phase begins. The magnetic energy is trans- 
ported through the inertial range up to the forcing scale 
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of the turbulence on roughly the local eddy timescale. 
From this scale, which is about the Jeans scale in the case 
of a star forming region, the field can be transported to 
even larger scales by outflows of stars or supernovae. 
How the small-scale dynamo operates in detail depends 
on the magnetic Prandtl number Pm, which is the ratio 
between kinematic viscosity v and magnetic diffusivity 
r\. With the hydrodynamic and magnetic Reynolds num- 
bers Re = VL/v and Rm = VL/r], where V is the typical 
velocity at the largest scale of the inertial range L, one 
can define Pm = Rm/Re. While the small-scale dynamo 
analytically is well studied in the limit of infinite Pm [Bi- 
ll], there are only a few studies for the case of Pm — > 
[9TTT2"]. Up to now, simulations have been restricted to 
the regime of 0.1 < Pm < 10 [13]. 

Nature features a broad range of magnetic Prandtl num- 
bers, reaching from about 10 12 in the primordial and 
present interstellar and intergalactic medium |14j to 
10~ 7 — 10~ 2 in the interior of planets and stars [T5I I16| . 
The magnetic activity in the Sun is well explained by a 
large-scale dynamo model [171119] . However, there are 
observational indications that also the small-scale dy- 
namo could play a role [20] [H] . In addition, the small- 
scale dynamo could provide a source for large-scale dy- 
namo action, especially if rotation is present [22] . More- 
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over, the regime of low magnetic Prandtl numbers is 
important in the case of liquid metal laboratory experi- 
ments, which also show the existence of a turbulent dy- 
namo (see, e.g., Nornberg et al. [13]). These provide 
an additional comparison for theoretical results from the 
Kazantsev theory or simulations. 

The efficiency of the small-scale dynamo depends on the 
type of turbulence. Most previous studies analyzed the 
case of ideal Kolmogorov turbulence [21], i. e. purely 
solenoidal turbulence. But astrophysical plasmas are 
usually highly compressible. For example, the energy and 
momentum input by supernova explosions lead to highly 
supersonic motions in the interstellar medium [2"51 - f2"T] . A 
similar interference is made for the accretion flow onto 
galactic disks or the convergent flows induced by spiral 
density waves 28 30J . Observations of the turbulent ve- 
locity spectrum confirm this argument [31) . Here we take 
into account the effects of different types of turbulence 
[31H35j ranging from incompressible Kolmogorov turbu- 
lence to highly compressible Burgers turbulence [55] . 
In this paper we describe phenomenologically how the 
small-scale dynamo operates at Pm — ¥ as motivated by 
the common stretch-twist-fold toy model, which was sug- 
gested for Pm — > oo. We summarize the concepts and the 
main equations of the Kazantsev theory, which analyti- 
cally describes the small-scale dynamo in the kinematic 
limit. As input, one requires the correlation function of 
the turbulent velocity field (e. g. see the model presented 
in Schober et al. [8]). We employ the WKB approxima- 
tion to solve the Kazantsev equation and find the growth 
rate of the magnetic field in the limit of small magnetic 
Prandtl numbers. The critical magnetic Reynolds num- 
ber Rm cr it, which needs to be exceeded for small-scale 
dynamo action, is the same for all magnetic Prandtl num- 
bers analyzed. We discuss the influence of Rm cr i t in the 
limit of small Pm. In the final section, we prove that the 
WKB approximation is valid for small to moderate Pm. 



II. PHENOMENOLOGY OF SMALL-SCALE 
DYNAMO GROWTH 

The small-scale dynamo converts kinetic energy from 
turbulent motions into magnetic energy. An illustrative 
model describing this process is the stretch-twist-fold dy- 
namo [37j . The stretching of a closed magnetic flux rope 
leads to amplification of the magnetic field strength, as 
the magnetic flux is a conserved quantity. Afterwards the 
rope is stretched, twisted and folded such that the orig- 
inal shape is regained. The shorter the turnover time of 
the turbulent eddies is, the faster the stretch-twist-fold 
mechanism proceeds and thus the faster the magnetic 
field is amplified. Intuitively the turnover time decreases 
with decreasing eddy length. 

In the limit of high magnetic Prandtl numbers the am- 
plification rate of the dynamo is most efficient on the 
smallest scale of the inertial range, i. e. the viscous scale 
l v = Re -1 ' L, where again $ is the slope of the tur- 



bulent velocity spectrum. 

During the transition from large to small magnetic 
Prandtl numbers, the resistive scale £ v = Rm _1 '^ +1 ' L 
becomes larger than the viscous one. The amplification 
then takes place at roughly £„, which lies in the iner- 
tial range of the turbulent velocity spectrum. Due to 
larger time scales of the turbulent eddies in the inertial 
range, we expect the small-scale dynamo to be less ef- 
ficient at low magnetic Prandtl numbers. While in the 
large Prandtl regime the hydrodynamical Reynolds num- 
ber regulates the dynamo, here the magnetic Reynolds 
number is the relevant quantity. 



III. ANALYTICAL DESCRIPTION OF THE 
SMALL-SCALE DYNAMO 

A. The turbulent velocity field 

A theoretical description of turbulence starts with the 
decomposition of the velocity field v into a mean field (v) 
and a turbulent component 6v: 

v = (v> + Sv. (1) 

The correlation of two turbulent velocity components at 
the positions ri and r 2 at the times t and s for a Gaussian 
random velocity field with zero mean, which is isotropic, 
homogeneous, and 5-correlated in time, is 

(6v i (r 1 ,t)6v j (r 2 ,s))=T ij {r)S(t-s) (2) 

with the two-point correlation function Ty(r) and r = 
\r± — r2 1 . The delta-correlation in time is a simplifying 
assumption, and its consequences should be explored in 
future studies. Following Batchelor [38 , Tjj(r) can be 
divided into a transverse part Tn and a longitudinal part 
Tl in the following way: 

T a( r ) = (*« - 7?) T n(0 + 7? r L(r). (3) 

We neglect here the effect of helicity, which would appear 
as an additional term in . Any turbulent flow can gen- 
erally be described by the relation between the velocity 
v{t) and the size I of a velocity fluctuation, 

v(£) cx P. (4) 

The power-law index i? varies from its minimal value of 
d = 1/3 for Kolmogorov theory [24 , i. e. incompressible 
turbulence, to Burgers turbulence [36], i. e. highly com- 
pressible turbulence, where i? gets its maximal value of 
1/2 [35]. 

We use the model for the correlation function of the tur- 
bulent velocity field from Schober et al. [8]. The longi- 
tudinal correlation function in the inertial range is moti- 
vated from the turbulent diffusion coefficient. We ensure 
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a continuous extension into the viscous range via an ap- 
propriate normalization. This leads to 



H(r) 




Re (i-tf)/(i+tf) (r) 



< r < 4 
4 < r < L 
L <r, 



(5) 

where 4 — L Re denotes the cutoff scale of the 

turbulence, i. c. the viscous scale, and L is the length 
of the largest eddies. In our model the transverse cor- 
relation function for the general slope of the turbulent 
velocity spectrum is 



VL 
3 



1 - i(#)Re 



(1-i5)/(1+i9) (r 







0+1 



< r < 4 
4 < r < L 
L < r, 

(6) 

with = (21 - 38i?)/5. The functional form of T N is 
based on the relation between the transversal and lon- 
gitudinal correlation function in the extreme cases of 
divergence-free (Kolmogorov) and rotation-free (Burg- 
ers) turbulence. 

B. Kazantsev Theory 

Like the velocity field, the magnetic field can be sepa- 
rated into a mean field (B) and a fluctuation part SB: 



B = (B) + SB. 



(7) 



Assuming that the fluctuating component <5B is a homo- 
geneous, isotropic Gaussian random field with zero mean 
like the velocity field, we can write down the correlation 
function as 



with the two-point correlation function 
M„(r,i) = (c5, 



(8) 



nr 



M N (r,t) + T -^M L (r,t). (9) 



As the magnetic field is always divergence-free the trans- 
verse and the longitudinal correlation function are related 
by 



Mk 



1_ d_ 

2r dr 



(r 2 M L ) 



where we have used that (r^j /r 2 )Mij 

(n/r^Mij = M N . 

The time derivative of M,, is 



(10) 
Ml and 



dM K 



d_ 
Of 



(SBiSBj) 



OB 



B 9B j 



d_ 
dt 



m) (b 3 )). 
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In the upper equation we can substitute the induction 
equation 

dB 



VxvxB-rjVxVxB, 



(12) 



where i] = c 2 /(Ana) is the magnetic diffusivity with the 
speed of light c and the electrical conductivity a, and the 
evolution equation of the magnetic mean field 



d(B) 
dt 



V x (v) x (B) - ?/ cff V x V x (B) 



(13) 



with the effective parameter rj c g = r\ + 2l(0). After a 
lengthy derivation [TO] this leads to 



dt 



2K dm M'{ + 2 



4«diff 



Miff 



Mr 



r \ r 



Tl 

r 



with 



Kdiff(r) =r? + T L (0)-T L (r). 



(14) 



(15) 



The prime denotes differentiation with respect to r. The 
diffusion of the magnetic correlations, /Miff j contains in 
addition to the magnetic diffusivity rj the scale-dependent 
turbulent diffusion T L (0) - T L (r). 

With the solution of Eq. ( 14 ) we can calculate Mn by us- 



ing the relation (10). Thus, we find the total correlation 
function of the magnetic field fluctuations Mjj, which is 
proportional to the energy density of the fluctuating part 
of the magnetic field, SB 2 /(8n). 

To separate the time from the spatial coordinates we use 
the ansatz 



M h (r,t) = 



1 



ip(r)e 



217 



Substitution of this ansatz in Eq. ( 14 1 gives us 



d 2 7 



+ U(r)ip(r) = -Tip(r). 



(16) 



(17) 



This is the Kazantsev equation, which is formally similar 
to the quantum-mechanical Schrodinger equation with a 
"mass" h 2 /(2n dm ) and the "potential" [40] 



Mr) 



Miff 



(n' dis ) 2 | 2 Kdi ff | 27* | 2(T L - T N ) 
4Kdiff t 2 r r 2 



(18) 



C. Formal Solution of the Kazantsev Equation in 
the WKB Approximation 

For the solution of the Kazantsev equation we use the 
WKB approximation. To use the standard formulation of 
this method, we have to make some substitutions. Def- 
inition of a new radial coordinate x with r = e x leads 
to 
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Next we eliminate the first-derivative terms through the 
substitution 



ijj{x) = e x/2 9(x), 



to obtain 



d 2 9(x) 



dx 2 



+ p(x)9(x) = 



with the definition 



D(x)= ir + u(x)}e 2 * 1 



(20) 



(21) 



(22) 



The WKB solutions of Eq. (21) are linear combinations 
of 



) = ^Tjcxp (±i jf y/p^dx'^ , (23) 



0(x 



where x\ is the first root of the p function p(x). The 
second derivative of ( 23 1 is 



'(x) 



1 



Ap 2 16 p 3 



where now the prime denotes d/dx. This equation results 
in the Kazantsev equation (17) if 



1/0*01 « i. 



with 



f s (j/y 

Ap 2 16 p 3 



(25) 



(26) 



From the shape of the p-function we conclude that the 
solutions between the two roots of p(x), i. e. X\ < x < X2, 
are oscillatory. 

The condition for the eigenvalues T is [4Tj 



yjp(x')dx' = 1 7r 



(27) 



for different excitation levels n € N. In this work we 
concentrate on the lowest mode n = 0, which has the 
largest growth rate. 




FIG. 1. (Color online) Dependence of the potential on the 
dimensionless parameter y = r/L for Kolmogorov ($ = 1/3) 
and Burgers ($ = 1/2) turbulence at Re = 10 s . We choose 
different magnetic Reynolds numbers Rm = 10 4 , Rm = 10 5 
and Rm = 10 6 , resulting in the Prandtl numbers Pm = 10~ 4 , 
Pm = 10~ 3 and Pm = 10" 2 . The viscous scale l v depends 
on the type of turbulence and the Reynolds number. For 
Kolmogorov turbulence l v = Re _3,/4 L; for Burgers turbu- 
lence i v = Re _2 ^ 3 L. The resistive scale is l v — Rm^ 3 ' 4 ! for 
Kolmogorov and £ v = Rm" 2/3 L for Bur gers turbulence. A 
magnetic Reynolds number 10 x is indicated in the resistive 
scale as (Re appears only in the viscous range.) 



IV. GROWTH RATE IN THE LIMIT OF SMALL 
MAGNETIC PRANDTL NUMBERS 

We are interested in bound eigenfunctions of the 



Kazantsev equation (21), which have corresponding real 
eigenvalues, i. e. growth rates. For this we require part 



of the potential ( 18 ) to be negative. 
In Fig. [I] we show the normalized potential ( 18 ) as a func- 
tion of y = r/L for Kolmogorov and Burgers turbulence. 
We choose a Reynolds number of 10 8 , which is a typi- 
cal value for example for the interior of planets |16j and 
primordial halos [14] . The different lines correspond to 



different magnetic Reynolds numbers of 10 4 , 10 5 and 10 6 
and hence represent magnetic Prandtl numbers of 10~ 4 , 
10~ 3 and 1CP 2 , respectively. The crucial discrepancy to 
the contrary limit of large Pm is that the potential only 
has a negative part in the inertial range (i. e. the range 
between i v and L indicated in the figure). Thus, there 
are only real positive eigenvalues of the Kazantsev equa- 
tion (17) in this range. 



With our model for the correlation function of the tur- 
bulent velocity field, Eqs. ^ and ([6]), thep function (22) 
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in the inertial range is 

p(y) 



^135 + Rmy (60 yT 



20(3 + Rmy 1 +»Y 
-o(0) Rmy 1+2€ 

+2/ (25-6(tf) + 10Rm?/ 2 f))) 



where 



r 



L 
V 



(28) 



(29) 



is the normalized growth rate, and we use the abbrevia- 
tions 



a(tf) = tf(56 - 103i?) 
b(tf) = ^(79 - 157tf). 



(30) 
(31) 



For the analytical determination of the zeros of p(y) we 
use the approximations 



Pi(y) 



3 Rm y 1+ o (a(i?) Rm y 1+i) + 2 b{$) - 50) - 405 



20 (3 + Rm?/ 1 +' ? ) z 
which is valid for r — > 0, and 

P2(y) = ^ Q — 3y r, 



(32) 



(33) 



where we leave out the constant terms in ([28]). We show 
p(y) as well as the two approximations in Fig. [2] for the 
exemplary case of Re = 10 8 and Rm = 10 5 . 
By using pi (y) we find for the first zero of p(y) approx- 
imately 



Ui = 



Rm 



(34) 



where we defined 



25 + Jl35a0&) + (6(t?) - 25) 2 - b(d) 
^ = ^ W) • ^ 



With P2(y) we find the second zero 



2/2 



20 r 



(36) 



The eigenvalue can be determined approximately by the 
equation 



L v 2 ' 



(37) 



Here we use the approximative function p2 (y) instead of 
the full function p(y) in order to find an analytical solu- 
tion of the integral. Note that the scaling of the abscissa 
in Fig. [5] is logarithmic and thus p^ (y) is a good approxi- 
mation of p(y) for y > y±. The value of the integral does 
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FIG. 2. (Color online) The p function (red curve) depending 
on the normalized scale parameter y — r/l in the viscous, 
inertial and large-scale range. We indicate the viscous scale 
l v and the forcing scale L as vertical lines. For the plot, we 
use our result for the growth rate in the limit of small Pm, 
Re = 10 8 and Rm = 10 5 . The dashed green line reefers to 
the approximation function pi and the dotted purple line to 
P2- The upper panel shows p(y) for Kolmogorov turbulence, 
the lower panel for Burgers turbulence. 



not change by much due to this simplification. 

We can solve the resulting equation from (37) with the 

ansatz 



r = a Rm i+* . 



(38) 



This is motivated by the result of Schober et al. [S] 
in the limit of large magnetic Prandtl numbers: F cx 
R e 0--&)/(i+v) , Here the amplification process takes place 
at the viscous scale, which depends on the hydrody- 
namical Reynolds number. As mentioned above in the 
limit of low magnetic Prandtl numbers the dynamo 
operates mainly on the resistive scale, which depends 
on Rm. Thus, our ansatz is to replace the hydrody- 
namic Reynolds number by the magnetic one (see also 
e.g. Boldyrev and Cattaneo [4"2"]). 
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With (|38|) we find for the solution of (|37|) 

' ' " ' ' s(i?) - 20c(^)t+* a+ 



0-1 V 5 

log (4 75 \/ RmT 



log (i 



c(0) 
Rm 



+ \Ja{$) - 20c(tf)^+l 



7T 

2' 



(39) 



As we assume the pre-factor of the growth rate to be very 
small, i. e. a < 1, we use a(<&) > 20 c^ 1 -^ 1 ^ a to 
approximate (39) as 



Rm 



log (4 y/5 V 

Rm CV 



The solution of this equation can easily be found: 



a(tf) 



c (^)i + i> exp 



3o(i?) 



7T (l? - 1) 



7T 
2' 

(40) 



(41) 



We list results for the normalized growth rate of the 
small-scale dynamo in the limit of low magnetic Prandtl 
numbers for exemplary types of turbulence in Tab.[TJ For 
comparison we also list the results for Pm — > 00 from 
Schober et al. [8]. Moreover, we present the normalized 
growth rate in Fig. [3| The lower abscissa shows the de- 
pendency on the magnetic Reynolds number, which is 
valid for any hydrodynamic Reynolds number Re 3> Rm. 
The upper abscissa in Fig. [3] shows the dependency on 
the magnetic Prandtl number for a fixed Re of 10 20 . We 
present the results for different types of turbulence re- 
ported in the astrophysical literature [2 



V. 



CRITICAL MAGNETIC REYNOLDS 
NUMBER 



For the onset of the small-scale dynamo the magnetic 
Reynolds number needs to exceed a critical value Rm cr j t . 
We determine the latter by setting the growth rate in 
(22) equal to zero and solving Eq. (27) for Rm. As the 



p-function in the inertial range only depends on Rm, but 
not on Re, it is independent of the magnetic Prandtl 
number. We list the numerical results for Rm cr ; t in Tab. 
|TJ The critical magnetic Reynolds number increases with 
increasing compressibility. 

In the limit of large Pm the critical magnetic Reynolds 
number is not necessarily the dominant restriction, as 
Re > 10 3 is required for turbulent flows. As Rm ^> Re 
for large Pm, Rm needs to be much larger than 10 3 , which 
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FIG. 3. (Color online) The normalized growth rate f as a 
function of the magnetic Reynolds number Rm (lower x axes) 
and the magnetic Prandtl number Pm (upper x axes). The 
results shown for the lower abscissa are only valid for small 
Pm, i. e., Rm -C Re, while we used a fixed Re of 10 s for 
the upper abscissa. We present different types for turbulence, 
indicated by the slope of the turbulent velocity spectrum 
K41 [24], SL94 [32], BNP02 [33], L81 [31], FRKSM10 [34] 
(sol: solenoidal forcing; comp: compressive forcing), OM02 
[35] and B48 M\- 



is larger than the critical magnetic Reynolds number. 
In the opposite limit of small Pm we have the case of 
Rm <C Re. For low hydrodynamic Reynolds numbers, 
the magnetic Reynolds number can fall below Rm cr it , and 
the small-scale dynamo can not operate. 
We note that, contrary to our results presented here, 
Iskakov et al. [13] found a weak dependence of the criti- 
cal magnetic Reynolds number on the magnetic Prandtl 
number. Therefore, it needs to be explored further 
whether the discrepancy in our results is due to approxi- 
mations in the Kazantsev model or if it is a result of the 
relatively narrow inertial range in numerical simulations. 
The latter provides a restriction on the number of turbu- 
lent eddies resolved in the box, and thus on the overall 
statistical sampling of the dynamics. In future studies, it 
would thus be desirable to explore this behavior at higher 
resolution in numerical simulations, and by relaxing the 
assumption of the Kazantsev theory in analytical studies. 



VI. COMPARISON WITH NUMERICAL 
SOLUTION 

Bovino, Schleicher, and Schober [44] solved the 



Kazantsev equation (17) numerically with the Numerov 
algorithm. They used the same form of the correlation 
functions of the turbulent velocity, i. e. Eqs. ^ and (|6j). 
A comparison of their result with our analytical solution 
is shown in Fig. [4] at a fixed Reynolds number of 10 14 . 
The comparison shows excellent agreement between 
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Model and reference 





T (Rm<Re) 


T (Rm>Re) 


Rnlcrit 


Kolmogorov |24j 


1/3 


0.027 Rm 1/2 


1.03 Re 1/2 


« 107 


Intermittency of Kolmogorov turbulence [32] 


0.35 


0.027 Rm 48 


0.94 Re ' 48 


w 118 


Driven supersonic MHD-turbulence |33| 


0.37 


0.026 Rm 46 


0.84 Re 46 


w 137 


Observation in molecular clouds 1311 

V^1_/.J^1 ill lilVJl^^ U1U,1 V^WLILIU ("^| 


0.38 


0.025 Rm ' 45 


0.79 Re ' 45 


« 149 


Solenoidal forcing of the turbulence [34] 


U.4o 


u.uiy nm 


u.04 ±ie 


~ 007 


Compressive forcing of the turbulence [34] 


0.47 


0.012 Rm 36 


0.34 Re 36 


w 697 


Observations in molecular clouds [35] 


0.47 


0.012 Rm ' 36 


0.34 Re ' 36 


w 697 


Burgers [36] 


1/2 


0.0054 Rm 1/3 


0.18 Re 1/3 


w 2718 



TABLE I. The normalized growth rate of the small-scale dynamo F in the limit of small magnetic Prandtl numbers (Rm<Re) . 
For comparison, we present also P for large magnetic Prandtl numbers (Rm3>Re). We show our results for different types of 
turbulence, which are characterized by the exponent $ of the slope of the turbulent velocity spectrum, v(£) oc £ . The extreme 
values of •& are 1/3 for Kolmogorov turbulence and 1/2 for Burgers turbulence. 




K41: = 1/3 
L81: i9 = 0.38 
FRKSMlO(sol): § = 0.43 
FRKSMlO(comp): # = 0.47 
B48: = 1/2 



Pm. This is probably caused by the approximations we 
made from Eqs. (39 1 and (40). 



10° 

Pm 



FIG. 4. (Color online) The normalized growth rate f as a 
function of the magnetic Prandtl number Pm. The results are 
shown for a fixed Reynolds number Re = 10 14 . We present 
the solutions from numerical integration of the Kazantsev 
equation by Bovino, Schleicher, and Schober [44] indicated as 
dashed lin es. The analytical solutions in the limits of small, 
i. e., Eq. (38 1, and large Prandtl numbers [H] are shown by 



the solid lines. We present different types for turbulence, indi- 
cated by the slope of the turbulent velocity spectrum K41 
[24], L81 [31], FRKSM10 [34] (sol: solenoidal forcing; comp: 
compressive forcing), and B48 36 . 



the numerical and the analytical solutions in the limit of 
small (f oc Rjn( 1_ *)/( 1+ *)) and large magnetic Prandtl 
numbers (f oc Re^-^/t 1 ^ |V We find that the range 
where our analytical solution ( [38] ) can be used is not re- 
stricted to Pm <C 1 but is also applicable in the regime 
Pm w 1 for all types of turbulence. We see a minor offset 
between our solutions and the numerical ones for small 



VII. CONCLUSION 

In this paper, we used the Kazantsev theory to de- 
termine the growth rate of the small-scale dynamo in 
the limit of low to moderate magnetic Prandtl num- 
bers. We found that the growth rate is proportional to 
Rm' 1_, '"' 1+ ''', where Rm is the magnetic Reynolds num- 
ber and $ is the slope of the turbulent velocity spectrum 
in the inertial range. The critical magnetic Reynolds 
number for small-scale dynamo action Rm cr it ranges from 
roughly 100 for Kolmogorov turbulence (z9 = 1/3) to 2700 
for Burgers turbulence ($ = 1/2). These values are the 
same for large and low magnetic Prandtl numbers. How- 
ever, for large Pm Rm cr i t provides no strong constraint, 
as here Rm 3> Re and Re > 10 3 for turbulence. In the 
limit of small Pm, where Rm <C Re, the critical magnetic 
Reynolds number is more important. We derived our re- 
sults employing the WKB approximation, which we have 
shown to accurately solve the Kazantsev equation in the 
limit of small magnetic Prandtl numbers (see Appendix 
|A| ) . Numerical integration of the Kazantsev equation pre- 
dicts a smooth transition of the growth rate between the 
small and large Prandtl number regime (see Fig. [4] and 
[44] ) . The analytical solutions for Pm <C 1 and Pm ^> 1 
cover a broad range of possible magnetic Prandtl num- 
bers and are also applicable at Pm w 1. This helps us to 
better understand numerical simulations, which, because 
of the limited resolution that can be achieved, are bound 
to magnetic Prandtl numbers around unity. 
We gained important results for the small-scale dynamo 
from the Kazantsev theory, which are summarized in Ta- 
ble [Tj However, one has to be careful with the indicated 
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numerical values as these show only the expected trends 
resulting from our assumptions. The growth rates as well 
as the critical magnetic Reynolds numbers can change 
if we take additional physical mechanisms into account, 
such as helicity or the finite correlation time of the turbu- 
lent eddies. Future highly resolved numerical simulations 
will provide a basis for comparison and help to deter- 
mine how strongly these additional effects may change 
the properties of the turbulent dynamo. 
With our calculations we show that small-scale dynamo 
action is not restricted to the regime of large magnetic 
Prandtl numbers, but it also occurs at low and moder- 
ate Pm. This can be used to analyze the evolution of 
the small-scale magnetic field in various physical envi- 
ronments in more detail. Rapid amplification of mag- 
netic seed fields is necessary to explain the strength of 
magnetic fields in the present-day Universe. 
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Appendix A: Validity of the WKB Approximation 

In this section we show that the WKB approximation, 
which we use for solving the Kazantsev equation, is valid 
in the limit of small magnetic Prandtl numbers. There- 
fore, we have to analyze (26) for the inertial range of the 



turbulence spectrum. As / is a function of the distance 
x, i. e. r, we have to evaluate it on the characteristic scale 
in which we are interested. With the main amplification 
occurring at the minimum of the potential, it is rational 
to use this scale. 



In Fig. [5] we show /(Pm) on the scale of the potential 
minimum for different Reynolds numbers, Re = 10 6 , 
Re = 10 8 and Re = 10 10 . We choose the example of 
Kolmogorov turbulence for the discussion. The test of 
validity is of course similar for other types of turbulence. 
The magnetic Reynolds number needs to exceed Rm cr i t : 



Rm = Pm Re > Rm c 



(Al) 



a 



0.15 
0.1 
0.05 


-0.05 
-0.1 
-0.15 




Re : 
Re : 
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10 s 



Re = 10 1 



1 i. 




-20 



-15 -10 

logio(Pm) 



FIG. 5. (Color online) Test of the validity of the WKB ap- 
proximation for Kolmogorov turbulence. The function ( 26 1 is 



shown depending on the magnetic Prandtl number Pm for dif- 
ferent hydrodynamic Reynolds numbers Re = 10 6 , Re = 10 s 
and Re = 10 10 . We evaluated /(Pm) at the minimum of 
The WKB approximation is valid for 



the potential ( 18 1 
/(Pm) — > 0. The critical magnetic Reynolds number gives 
a further restriction (see text) leading to possible values for 
Pm of Pm > 10" 4 (Re = 10 6 ), Pm > 10" 6 (Re = 10 s ) and 
Pm > 10 -8 (Re = 10 10 ), which we indicated by the corre- 
sponding arrows. 



For Kolmogorov turbulence Rm cr i t ~ 10 2 . Thus, for ex- 
ample in the curve with Re = 10 8 in Fig. [5] the threshold 
is only exceeded for Pm > 1CK 6 . In this regime the WKB 
approximation is perfectly valid. In Fig. [5] we use arrows 
to indicate the regimes, where the small-scale dynamo 
can operate. 

In principle, Fig.[5]states that our approximation is valid 
also for larger Pm up to Pm — > 00. However, at some 
point we face the problem in which the potential gets 
negative also below the viscous range. For Pm — > 00, 
the negative part of the potential in the viscous range 
clearly dominates. As we do not account for this range 
in the calculations above, our results are only valid for 
sufficiently small Pm. 
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